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ABSTRACT 



1RXS J101015. 9-31 1909 is a galaxy located at a redshift of z = 0.14 hosting an active nucleus (called AGN) belonging to the class of bright 
BL Lac objects. Observations at high (HE, E > 100 Me V) and very high (VHE, E > 100 GeV) energies provide insights into the origin of 
very energetic particles present in such sources and the radiation processes at work. We report on results from VHE observations performed 
^ ■ between 2006 and 2010 with the H.E.S.S. instrument, an array of four imaging atmospheric Cherenkov telescopes. H.E.S.S. data have been 
5-H ■ analysed with enhanced analysis methods, making the detection of faint sources more significant. VHE emission at a position coincident with 
1RXS J101015. 9-31 1909 is detected with H.E.S.S. for the first time. In a total good-quality livetime of about 49 hours, we measure 263 excess 
counts, corresponding to a significance of 7.1 standard deviations. The photon spectrum above 0.2 TeV can be described by a power-law with a 
photon index of T = 3.08 ± 0.42 stat ± 0.20 sys . The integral flux above 0.2 TeV is about 0.8% of the flux of the Crab nebula and shows no significant 
variability over the time reported. In addition, public FermifLAT data are analysed to search for high energy emission from the source. The 
Fermi/LAT HE emission in the 100 MeV to 200 GeV energy range is significant at 8.3 standard deviations in the chosen 25-month dataset. UV and 
X-ray contemporaneous observations with the Swift satellite in May 2007 are also reported, together with optical observations performed with the 
ATOM telescope located at the H.E.S.S. site. Swift observations reveal an absorbed X-ray flux of F ( o.3-7)kcV = 1.04^^ x 10~" ergcrrT 2 s _I in the 
0.3 - 7 keV range. Finally, all the available data are used to study the multi-wavelength properties of the source. The spectral energy distribution 
(SED) can be reproduced using a simple one-zone Synchrotron Self Compton (SSC) model with emission from a region with a Doppler factor of 
30 and a magnetic field between 0.025 and 0. 16 G. These parameters are similar to those obtained for other sources of this type. 

Key words, y-rays - Galaxies: active - BL Lacertae : individual : 1RXS J101015.9-31 1909 
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1. Introduction 




BL Lac objects are characterised by rapid variability in all en- 
ergy ranges, and often display jets with apparent superluminal 
motions. Their extreme properties are thought to be related to 
the relativistic bulk motion of the emitting region at small an- 
gles to the line of sight of the observer. In addition, these ob- 
jects show highly polarized emission and no or only weak emis- 
sion lines. The observed broadband Spectral Energy Distribution 
(SED) of BL Lacs is often comprised of two bumps, one peak- 
ing at lower (radio to X-ray), the other peaking at higher (above 
X-ray) energies. In leptonic scenarios, the lower energy compo- 
nent is generated by synchrotron emission of relativistic elec- 
trons moving inside the jet. The higher energy component is due 
to the inverse Compton scattering of the electrons off the pho- 
tons of the self-gene rated synchrotron photon field (SSC mod- 
els, see for instance Marsch er & Gear I I1985I) . or off the pho- 
tons provided externally by other region s of the source (External 
Comp ton or EC models, see for instance lbermer & Schlickeiser I 
1 1993b . The VHE y-ray emission in hadronic scenarios can also 
be explained b y the interactions of relativistic proton s with am- 
bient photons (iMannheiml l 1 993b or magnetic fields dAh aronian 
2000). Depending on the position of the synchrotron component, 
BL Lacs are subdivided into Low-frequency peaked (LBL) if 
the maximum of the emission is in the infrared band, and High- 
frequency peaked (HBL) if the emission is peaked in the UV/X- 
ray band. 

1RXSJ101015.9-311909 belongs to the ROSAT All Sky 
Survey Bright Source Catalog (RA SS/BSC) of soft (0.1-2 keV) 
X-ray sources (Voges et al. 1119991) . with a flux of 2.9 x 10" 11 
erg cm" 2 s" 1 . It is located at a position of (arj20oo> feooo) = 
(10 h 10 m 15.9 s ,-31°19'09 "i!l and has a redshift of z = 0.14 (see 
Piranomonte et al. 2007 for both measurements). The source 
is present in the NRAO VL A Sky Survey (NVS S) catalogue 
of radio sources at 1.4 GHz dCondon et al. Ifl998h . which lists 
its flux density as 73.5 + 2.7 mJy. A radio flux of 89.5 + 
3.6 mJy in the 843 M Hz band has be en measured by the 
SUMSS radio survey dBock et al. 1 1201 lh . Due to its extreme 
value of the X-ray to radio flux ratio and its high X-ray flux, 
1RXSJ101015.9-3 11909 passed the criteria fo r inclusion in 
the Se dentary Multi-Frequency Survey catalogue (Giommi et al.l 
l200l . This catalogue specifically selected HBLs and thus pre- 
sented an obvious choice for the extension of the list of VHE 
BL Lac c andidates 1 RXS J 1 1 1 5 .9 -3 1 1 909 also fulfilled the 
criteria proposed in ICostamante & Ghisellinil (12002). where BL 
Lac candidates are considered interesting targets if they exhibit 
high levels of both X-ray and radio emission. 

Following these indications, observations of this source with 
H.E.S.S. started at the end of 2006, yielding the discovery of 
y-ray emission from 1RXS J101015.9-31 1909 (see Sec. re- 
ported here. By combining this information with other multi- 
wavelength data, the properties of the detected emission and 
its physical implications are discussed. The HE emission of the 
source has been studied with Fermi /LAT public data between 
100 MeV and 200 GeV and results are reported here in Sec. 13.11 
Analysis of data at lower energy bands is carried out to under- 
stand the emission from this source: Swift data (from the XRT 
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Fig.l. Image of 1RXS J101015.9-31 1909 in right ascension 
and declination (J2000) of the y-ray excess found by H.E.S.S. 
oversampled with the 68% containment radius of the point 
spread function (0.11° for these analysis cuts). The cross repre- 
sents the nominal position of the source. The inset on the lower 
left shows the expected y-ray excess distribution from a point- 
like source. 

and UVOT telescopes) are analysed and discussed in Sec. 13.21 
and !3.3l and optical data from the ATOM ( Automatic Telescope 
for Optical Monitoring. lHauser et al. 2004) telescope located on 
the H.E.S.S. site and taken mostly contemporaneously to the 
H.E.S.S. data, are analysed and discussed in Sec. 13.41 Finally, 
in Sec. |4] all the available data are used to study the global SED 
of the source in the context of a simple SSC scenario. 

2. H.E.S.S. observations and results 

H.E.S.S. is an array of four imaging Cherenkov telescopes lo- 
cated in t he southern hemispher e in the Khomas Highland of 
Namibia ( Aharonia n et al. 2006), that detects cosmic y-rays in 
the 100 GeV to 100 TeV energy range. Each of the telescopes is 
equipped with a segmented mirror of 107 m 2 area and a cam- 
era composed of 960 photomultipliers covering a large field-of- 
view (FoV) of 5° diameter. The stereoscopic system works in a 
coincidence mode, requiring at least two of the four telescopes 
to trigger the detection of an extended air shower. The trigger 
threshold, defined as the p eak of the di fferential y-ray rate for a 
Crab-like source at Zenith dFunk et al. 1120041) . is about 100 GeV 
and increases with increasing zenith angle. 

Observations of 1RXS J101015.9-31 1909 were carried out 
with H.E.S.S. in a campaign of 64 hours of observation time 
between 2006 and 2010. These cover a range of zenith angles 
between 8° and 28°, giving an average zenith angle of 12.9°, 
with a pointing offset of 0.5° relative to the nominal position of 
the source (see Tab.Q]for all details). The data from a total high- 
quality livetime of ~ 48.7 hours (after hardware and weather 
quality selection cri teria were applied with a procedure similar 
to that described in [Aharonian et al."1 l2006) have been analysed 
to search for emission at the nominal position of the source. 

The analysis of the > 100 GeV y-ray emission from this 
AGN is carried out with the analysis procedure described in 
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Table 1. Summai-y of good-quality data of H.E.S.S. observations of 1RXS J101015.9-31 1909 over the years 2006 - 2010. 



year 


MJD (start) 


MJD (end) 


N. 

1 'runs 


LT 


zen 


A'on 


Nopp 


N y 


cr 


2006 


54090.09 


54090.11 


1 


0.43 


12.7 


17 


126 


5.5 


1.5 


2007 


54142.95 


54238.79 


35 


14.52 


11.7 


551 


4835 


111.5 


4.9 


2008 


54475.07 


54535.92 


12 


5.37 


10.1 


136 


1291 


18.6 


1.6 


2009 


54832.06 


54976.79 


36 


15.62 


14.4 


457 


3968 


96.2 


4.6 


2010 


55265.90 


55299.84 


29 


12.80 


13.6 


255 


2466 


30.8 


1.9 


Tot. 


54090.09 


55299.84 


113 


48.70 


12.9 


1416 


12686 


262.7 


7.1 



Notes. The columns represent the year in which the source has been observed, the start and end date of observations in MJD, the number of 
good-quality runs available N mm , the corresponding exposure time in hours (LT), the mean observation zenith angle zen in degrees, the number of 
ON- (A'on) and OFF-source (A'off) events, the number of excess events N y and the significance of the detection in units of standard deviations a. 
The observation offset and the background normalization factor a (see text) for all datasets presented in the table are 0.5° and 0.09, respectively. 



iBecherini et al. I (1201 ll) . where an enhanced lo w-energy sensitiv- 
ity w ith respect to standard analysis methods (lAharonian et al. I 
2006) is achieved. This new analysis method is based on a 
multivariate signal-to-background discrimination procedure us- 
ing both previously-known and newly-derived discriminant vari- 
ables which depend on the physical shower properties, as well 
as its multiple images. In order to have a lower threshold for this 
source, the analysis configuration with a charge value of 40 pho- 
toelectrons has been used as a minimal required total amplitude 
for the cleaned and parametrized image in each telescope. 

The VHE y-ray emission from the BL Lac object 
1RXS J101015.9-311909 is detected using the Re flected back- 
ground modelling method (lAharonian et al. 1120061) with a statis- 
tical significance! of 7.1 standard deviations. The significance 
of the detection is represented by an excess of 263 counts at 
the nominal position of the source, the total number of ON- and 
OFF-source events being Non = 1416 and Noff = 12686, re- 
spectively, with a background normalization factoiQa' = 0.09. 

The VHE y-ray excess image obtained with th e Ring back- 
ground modelling method (lAharonian et al. |[2006l) is shown in 
Fig.[TJ while Fig.|2]shows the ON-source and normalized OFF- 
source angular distributions (9 2 ) for all H.E.S.S. observations: 
the background is rather flat, as expected at very small 2 , and 
there is a clear excess at small values of 6 2 , corresponding to the 
observed signal. A fit to the excess events of a point-like source 
model convolved with the H.E.S.S. point-spread-function (PSF) 
yields a position a J20 oo = 10 h 10 m 15.03 s + 3.77^ ± 1.56| ys and 
feooo = -31°18'18.4" ± 41.6'/^ + 20^, consistent with the po- 
sition of the radio and X-ray source (see Fig. |3}- The 3 <x upper 
limit to the intrinsic source extension calculated at the best fit 
position is 3.4'. 

The time-averaged differential VHE y-ray spectrum of 
the sourc e, derived using th e forward-folding technique de- 
scribed in [Piron et al. I (|2001), is presented in Fig. |4] The spec- 
trum is well fitted by a power-law function dN/dE = <po X 
(E/lTeVy r with a normalization of </>o = (1.87 ± 0.66 sta t ± 
0.37 sys ) x lO-^crrrV^eV-'and photon index F = 3.08 ± 
0.42 stat + 0.20 sys . The differential flux at the decorrelation en- 
ergy (£dec = 0.51 TeV) is fe cc = (1.47 ± 0.31 stat + 0.29 sys ) x 
10~ 12 cnr 2 s -1 TeV _1 . The integral flux above the analysis thresh- 
old £ th = 0.2 TeV is (f>(E > E th ) = (2.35 + 0.64 stat + 0.47 sys ) x 

2 Calculated following Eq. (17) of lLi&Mal fp983l) . 

3 In the Reflected background method a is just the reciprocal of the 
number of OFF-source regions considered. 

4 For this analysis the threshold energy is defined as the energy at 
which the effective detection surface exceeds two hectares and where 
the energy bias is less than twice the energy resolution. 
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Fig. 2. The distribution of 9 2 for ON-source events and normal- 
ized OFF-source events centered at the nominal position from 
H.E.S.S. observations of 1RXS J101015.9-31 1909. In the inset: 
8 2 distribution of the excess at the fitted position of the AGN; the 
superposed line is a fit of the PSF to the data. The containment 
radius at 68% of the PSF for the given observation conditions is 
0.1 16° and is shown by the dashed vertical line. 



10 -12 cm~ 2 s _1 , corresponding to ~ 0.8% of the flux of the Crab 
nebula above the same threshold. No significant variability is 
detected; the integral flux is seen to be constant within errors 
over the H.E.S.S. dataset, as shown in Fig. [5] A fit of the period- 
by-periocQ light curve with a constant value yields a^- 2 /dof = 
11.39/9, with a probability of 25%. The measured normalized 
excess variance of 0.44 + 0.71 on the same light curve yields a 
99% confidence level upper limit on the f ractional varianc e^ of 
< 151 %, as calculated using the method of iFeldman & Cousins I 
(119981) . No variability can be seen either in other time binnings 
tested (year-by-year or run-by-run). All analysis results have 
been cross-checked and confi rmed with an independent method 
(Ide Naurois & Rolland 1 12009). which gives consistent results. 



5 A H.E.S.S. observing period is the period between two full moons. 

6 See Vaughan et al.] i2003l) for definitions of normalized and frac- 
tional excess variance. 
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Table 2. Spectral properties for the analysis of Fermi/LAT data. 
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Fig. 3. Upper panel. The three contour lines around the H.E.S.S. 
fitted position of 1RXS J 101015.9-3 1 1909 correspond to the er- 
ror contours at 1 , 2 and 3 standard deviations on the position 
evaluation. The red dashed line represents the 1 <x significance 
contour of the Fermi/LAT detection of the source. Bottom panel. 
Zoom on the H.E.S.S. fitted position, represented by the cross 
and the error contours, compared to the X-ray (Swift/XRT and 
ROSAT) and radio (NVSS) positions. 



3. Multi-wavelength observations 

3. 1. Analysis of Fermi/LAT data 

1RXS J101015.9-311909 has been associated with the object 
2FGL J 1 009 .7—312 3 in the Fermi/LAT second source catalogue 
(Abdo et aT1 l201lh . A Fermi/LAT data analysis is performed 
on the publicly available data, spanning the time interval from 
2008-08-04 (MJD 54682) to 2 011-01-01 (MJD 555 62), using 
the binned likelihood method (lAtwood et al. 1 120091) from the 
Science Tools package V. v9r23pl, following the procedure rec- 
ommended by the Fermi/LAT collaboratiorQ 

The isotropic model iso_p7v6source is used to account for 
both the extragalactic diffuse emission and residual instrumental 
background, while the spatial template gal_2yearp7v6_v0 is 





r 


TS 


<KE > End 


^dec 


100 


2.09±0.15 stat 


68.29 


11.31 ±3.83 stat 


1929 


300 


1.92±0.15 stat 


62.73 


2.54 ± 0.69 stat 


3259 


500 


1.82 + 0.15 stat 


59.82 


1.36±0.34 slal 


4306 


1000 


1.71±0.16 stat 


55.93 


0.70±0.16 8tat 


5863 



see http://fermi.gsfc.nasa.gov/ssc/data/analysis 



Notes. The columns correspond to the energy threshold in MeV, Efh, the 
photon index, T, the test statistic (TS), the integral flux above threshold 
4>(E > E t h) in units of lO^ph cm 2 s~' and the decorrelation energy (£dec) 
in MeV. The Fermi/LAT systematic uncertainty on the spectral index is 
10% at 100 MeV, decrea sing to 5% at 560 M eV and increasing to 10% 
at lOGeV and above, see lAbdo et al. I ( 1201 ll) . 



used to account for the contribution from the Galactic diffuse 
emission. 

Since 1RXS J101015.9-31 1909 lies at a Galactic latitude of 
20.05°, the centre of the region of interest (Rol) is taken 5° away 
in the North-East direction from its nominal position in order to 
minimize the contribution from the Galactic diffuse emission. 

In the analysis presented here, source-class events are con- 
sidered in a circular Rol of 10° radius, and the P7V6_S0URCE 
instrumental response functions were used. In order to account 
for the potential contamination of events from sources outside 
the Rol due to the large PSF at low energies, all the neighbour- 
ing 2FGL objects are included in the model reconstruction of the 
source up to a radius of 15°. 

Using the gtlike tool and assuming a power-law sha pe for 
the source spectrum, the Test Statistic (TS. iMattox et al. I ll996) 
of the binned likelihood analysis is 68.3, corresponding approx- 
imately to a 8.3 cr detection in the 100 MeV - 200 GeV energy 
range. The corresponding photon index is F = 2.09 + 0.15 sta t 
and the highest energy photon from the direction of the source 
(i.e., within the 95% containment radius of the PSF at the given 
energy) has an energy of 76.6 GeV. Other, more complex spec- 
tral shapes like a log-parabola or a broken power-law do not re- 
sult in a significant improvement of the fit, and thus the power- 
law spectral shape is used in the remainder of this paper. The 
resulting spectral slope under these assumptions is consistent 
with the value found in the 2FGL catalogue, which gives Y = 
2.24 + 0.14 stat . 

However, there is evidence for a dependence of the photon 
index on the chosen energy threshold in the data analysis as sum- 
marized in Tab. [2] The spectrum of the source tends to harden 
with an increasing low-energy cut, which could be an indication 
of a curved spectrum. Future observations with Fermi/LAT may 
enable a significant detection of a possible curvature of the spec- 
trum compared to a pure power-law. 

To further check these results, a test was performed by mod- 
elling the Galactic diffuse emission with a power-law spectrum, 
instead of using a constant flux normalisation for this compo- 
nent, as is usually recommended by the Fermi/LAT team. Such 
an energy-dependent spectrum for this component would be an 
indication for a mis-modelled Galactic diffuse emission in the 
analysis, and could affect the hardening tendency as a func- 
tion of the energy threshold reported in Tab. |2] When using a 
threshold of 100 MeV, the latter test results in a photon index 
of F = 0.07 + 0.01 for the Galactic diffuse component, while 
the spectral results for the AGN remain fully compatible with 
those reported in Tab. [2] This slight energy-dependence of the 
spectrum of the Galactic model just reflects the fact that the 
mechanism responsible for the HE emission from the Galaxy is 
not yet perfectly understood, but does not strongly affect our re- 
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suits. While at each of the energy thresholds the count map of the 
Rol exhibits a visible gradient due to the Galactic diffuse emis- 
sion, no such gradient is present in the residual map after sub- 
traction of the Galactic and extragalactic models and the 2FGL 
sources (including 1RXS J101015.9-31 1909), being rather flat 
within the counting error. This shows that the normalization of 
the Galactic diffuse emission is under control and well-modelled 
in this analysis. 

In the following, we will adopt the results of Fermi/LAT data 
analysis using the two energy thresholds of 300 MeV and 1 GeV, 
see Fig. [4] 

The choice of a 300 MeV threshold is made in order to min- 
imize a possible contamination at low energies from neighbour- 
ing sources and from the Galactic diffuse emission. This choice 
takes into account the tendency of the spectrum to harden with 
increasing energy threshold, while not losing too many source 
photons due to this cut. We choose 1 GeV as a second thresh- 
old in order to study how the evaluation of the Fermi/LAT slope 
affects the modelling of the overall SED (see Fig. |7J. 

The Fermi/LAT binned spectral points shown in Fig. [7] are 
computed by running gtlike in five contiguous energy bins, us- 
ing the model parameters from the likelihood fit on the en- 
ergy range 1 GeV - 200 GeV, where the spectral index of 
1 RXS J 1 1 1 5 .9 - 3 1 1 909 was fixed to the best value of T = 1.71 
(see Tab.|2]i. An upper limit on the flux in a given energy bin was 
computed if TS < 9. The resulting fluxes for all analyses can be 
found in Tab. [2] 

The Fermi/LAT light curves, for the two chosen threshold 
energies, are shown in Fig. [5] where the data are presented in a 
6-month binning: given the low photon statistics, no significant 
variability is found in the 25 months of data. This was checked 
using other time binnings ranging from 90 to 180 days. 

The Fermi/LAT position of 1RXS J101015.9-31 1909 has 
been optimized using the tool gtfindsrc, and the best 
fit was found to be at the position (aj2ooo, £j200o) = 
(10 h 09 m 49.51\ -31°24'21.9") which is fully consistent with the 
position reported in the 2FGL catalogue (~ 3' away). The lcr 
contour presented in Fig. |3]was derived from the TS map com- 
puted on the Rol, using the best-fit position of the source. 

3.2. Swift/XRT 

The X-ray Telescope (XRT) dBurrows et al l 120051) on board 
the y-ray burst mission Swift dGehrels et al.l |2004|) observed 
1RXSJ101015.9-311909 three times during 2007-05-17 and 
2007-05-18 (see Tab.[3]for the total exposure time available with 
Swift/XRT). The first and third observations were performed in 
photon-counting (pc) mode, while the second observation was 
performed in windowed-timing (wt) mode. Cleaned event files 
have been reduced using HEASo/q V. 6.7. Source spectra and 
lightcurves have been extracted using XSelect, V. 2 . 4a, and the 
spectral fitting has been performed using XSpec, V. 12.5.1. 
Response matrices and ancillary response files have been pro- 
vided by the Swift/XRT instrument team. The source count-rate 
is equal to 0.4 counts s _1 for the three observations. The pres- 
ence of a pile-up effect in the data has been checked following 
the prescriptions of the Swift/XRT instrument teanfl leading to 
the conclusion that it does not affect the observations. As no sig- 
nificant variability has been observed, the two spectra obtained 
in the pc-mode have been summed using mathpha, V. 4.1.0, 
and fitted together with the second observation spectrum. Data 



8 http://heasarc.nasa.gov 

9 http://www.swift.ac.uk/pileupthread.shtml 
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Fig. 5. Upper panel. Light curve of H.E.S.S. observations. The 
mean flux above 0.2 TeV per observing period (between two full 
moons) is shown as a function of the time in MJD. Only sta- 
tistical errors are shown. Upper limits at 99% confidence level 
are calculated when no signal is found and in this case the cor- 
responding negative fluxes are shown in magenta. The solid line 
represents the fit of a constant to the H.E.S.S. data. Middle panel. 
Light curves of Fermi/LAT observations for the E± > 300 MeV 
and F t h > 1 GeV thresholds in a 6-month binning. The first two 
flux points of the F t h > 300 MeV light curve and the first three 
flux points of the F t h > 1 GeV light curve have a TS < 9. Bottom 
panel. Light curve of ATOM observations with R and B filters, 
corrected for Galactic extinction assuming E B -v = 0.104 (case 
A). 



below 0.3 keV have not been included in the analysis[3 while 
the last significant bin is at « 7 keV. The spectra have been re- 
binned using grppha, V. 3 . . 1, in order to have a minimum 
of 10 counts per bin. The Galactic column density A^h has been 
fixed a t 7.79 x 10 20 cm 2 , as evaluated by Dick ey & LockmarJ 
<fl990h . 

A fit performed using a simple power-law function with 
Galactic absorption gives T = 2.15 + 0.06 and normalization 
factor C lkeV = (3.0 + 0.1) x lO^keV-'cirrV 1 0t 2 /dof = 
172/141). The fit is significantly improved (F-test probability 
equal to 4 x 10~ 6 ) if a broken power-law is assumed, as shown 
in Tab. [4] case A, where the best fit parameters for the two pho- 
ton indices, break energy and normalization are presented. The 
absorbed flux in the 0.3 - 7keV energy band is found to be 
(1.04+°-<*) x lO^'ergcm-V 1 . 



10 http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt/ 
S WIFT-XRT-C ALDB -09.v 1 6.pdf 
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Energy (TeV) 



Fig. 4. Left panel. Time-averaged VHE spectrum measured from the direction of 1RXS J 101015.9-31 1909. The bow-tie represents 
lcr confidence level error band of the fitted spectrum using a power-law hypothesis. The lower panel shows the fit residuals, i.e. 
(N y - Ntheo) /Mheo, where N y and iV t heo are the detected and expected number of excess events, respectively. Right panel. Fermi/LAT 
and H.E.S.S. bow-ties. The two Fermi/LAT bow-ties represent the > 300 MeV (dotted line) and the £ t h > 1 GeV (dashed line) 
spectral results. All the bow-ties (H.E.S.S. and Fermi/LAT) result from a forward-folding spectral analysis technique. 



Table 3. Swift observations available for 
1RXSJ101015.9-311909. 





ID 


Mode 


Start 


Exposure (s) 


obs. 1 


00030940002 


pc 


2007-05-17 


1744 


obs. 2 


00030940003 


wt 


2007-05-18 


790 


obs. 3 


00030940004 


pc 


2007-05-18 


1981 


tot 








4515 



Notes. In the photon-counting (pc) mode the entire charge-coupled de- 
vice is read out, while in the windowed-timing (wt) mode only the cen- 
tral rows of the camera are read, increasing the time resolution of the 
instrument. 



Table 4. Parameters of the two hypotheses under consideration 
for the fit of the Swift/XRT data. 



Case 


r 


Ebreak [keV] 


ClkeV 


-^H.frcc 


X z /d.o.f. 


A 


r. = i.s!- 
r 2 = 2.5$ 


1 4+0.5 


q 9+0.2 
3 -0.2 




144/139 


B 


r = 2.5-; 




O Q+0.4 

o.y_ n4 


9!t 


147/140 



Notes. T is the fitted photon index, and £b, M k ls the energy at which the 
break in the spectrum occurs. Case A represents the broken power-law 
hypothesis considering only the absorption in the Galaxy, while case B 
represents the power-law fit taking into account the absorption in the. 
Galaxy plus a second absorber (iVufr*, expressed in units of 10 20 cnT 2 ) 
located at the redshift of the source. The normalization C\ k e v is given 
in units of 10~ 3 keV~'cnr 2 s~' 



The break observed in the X-ray spectrum (see Fig. |6) can 
either be intrinsic or external, i.e. due to an additional absorp- 
tion component in the AGN host galaxy. For an extensive dis- 
cussion about this topic see lPerlman et al. I (120051) . where based 
on the analysis of XMM-Newton spectra of 13 different BL Lac 
objects, a discussion of the intrinsic or external origin of the ob- 




1( Y-12L_, 1 1 , ill 1 1 1 1 1— I 

1U 0.3 0.4 0.5 1 2 3 4 5 6 

Energy (keV) 

Fig. 6. The two Swift/XRT spectra for case A (galactic absorp- 
tion only) and B (additional absorption), corrected for the re- 
spective absorptions, are shown in this figure by the red points 
and the black triangles, respectively. For simplicity, only the 
spectrum from the pc-mode observations, rebinned for plotting 
purposes, is presented. 



served spectral curvature is given, concluding that the first hy- 
pothesis would be preferred. The hypothesis of an external origin 
of the break has been tested by fitting the Swift/XRT data with 
a power-law emission function inclu ding absorption by Galactic 
material (fixed at the value given bv lDickev & Lockmarf l990) 
plus a second absorber located at the redshift of the host galaxy 
with adjustable column density. The best fit in this case (case 
B in Tab. |4]i is statistically equivalent to the broken power- 
law, the evaluation of the second absorber column density being 
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Nn,&ee = 9^ x 10 20 cm 2 . This second absorber is, however, 
poorly constrained compared to the Galactic one. 

It should be noted that, given the relatively low redshift of 
the source, the location of the absorber cannot be constrained. 
In particular, the same absorption effect could be obtained by 
multiplying by a factor of ~ 2 the contribution of the Galactic 
absorption in the direction of the source. However, such a high 
value of the Galactic column density is not consistent with the 
range of Nh measu red in a circle of 1 ° aroun d the nominal posi- 
tion of the source (iDickev & Lockmanll 1990b . 

The deabsorbed X-ray spectra of the source assuming either 
an intrinsic break of the spectrum (corrected only for Galactic 
absorption), or an external one (corrected for both absorbers), 
are shown in Fig. [6] 



3.3. Swift/UVOT 

The Swift satelli te carries an Ultra- Violet/Optical 

Telescope (UVOT) (Romingetal. 2005), which observed 
1RXSJ101015.9-311909 simultaneously with XRT Six dif- 
ferent filters are available: V and B in optical and U, UVW1, 
UVM2 and UVW2, in the ultra-violet, in order of increasing 
frequency. Counts have been extracted in a 5" radius of aper- 
ture, and magnitudes and fluxes have been evaluated using 
uvotmaghist, V. 1.1. The correction for Galacti c extinction 
has been done following iRoming et al.l (120091) . assuming 
Eb-v = 0.104 and 0.224 for case A and B respectively, where 
Eb-v is the difference of the total extinction in the B and V 
filters. The evaluation of Eg-v has been done using Nh = 7.79 
and 16.79 x 10 20 cm" 2 , for case A and B, respectively, and 
Nu/Eb -v = 7.5 x 10 21 cm" 2 , as given in iJenkins & Savage! 

No significant variability in the data is observed. 
Therefore, the mean flux values measured by UVOT for each 
filter are used for the study of the SED (see Sec. |4}. 



3.4. ATOM 

ATOM dHauser et alj|2004l) is a 75-cm optical telescope located 
at the H.E.S.S. site. 1RXS J101015.9-3 11909 has been regu- 
larly observed with ATOM from January 2008 to June 201 1. On 
2008-02-01 th e source has be en observed with B, R and I filters 
(as defined by Bessell 1990) while the rest of the observations 
have been performed in the B and R bands, only. The fluxes have 
been determined using a 4" radius of aperture. The evaluation of 
the errors includes the uncertainty of the absolute calibration. As 
shown in the lower panel of Fig. [5] the source varies significantly 
in both the B (79% +11%, evaluated as the difference between 
the highest and the lowest flux measured over the mean value) 
and R (62% + 8%) energy bands. The measured normalized ex- 
cess variances of (0.014 + 0.005) B and (0.011 + 0.002)« con- 
firm that there is a significant variability in both the blue and red 
bands, respectively. The lowest detected variability time-scale is 
roughly 1 day, corresponding to the minimum time between two 
different observations. The measured I band spectral flux den- 
sity, not shown in Fig. [5] is (7.84+0.33)x 10" 16 erg cm" 2 s" 1 A" 1 . 
For the SED (Sec.|4]i, the mean flux obtained from ATOM data, 
corrected for Galactic extinction (again using Eb-v = 0.104 and 
0.224 for case A and B, respectively), is considered; the error 
bars show the flux variability range observed. 



Table 5. Parameters used for the SSC modelling of the SED of 
1RXS J101015.9-31 1909 and derived physical quantities. 





Fermi >30 oMeV 


Fermi >1GeV 




Case A 


Case B 


Case A 


CaseB 


/ g, break 


1.08 x 10 5 


9.0 x 10 4 


1.01 x 10' 


7.8 x 10 4 


a\ 


2.2 


2.2 


2.0 


2.0 


K 


5.67 x 10 4 


3.37 x 10 2 


6.6 x 10 3 


4.5 x 10 1 


U e 


5.37 x 10" 2 


3.15 x 10" 4 


3.42 x 10" 2 


2.24 x 10" 4 


B 


0.16 


0.025 


0.16 


0.025 


u B 


1.02 x 10" 3 


2.49 x 10" 5 


1.02 x 10" 3 


2.49 x 10" 5 


uJu B 


52.7 


12.7 


33.6 


9.0 


R 


2.37 x 10 15 


7 x 10 16 


2.37 x 10 15 


7 x 10 16 


^"var 


0.8 


24.7 


0.8 


24.7 




7.6 x 10 42 


4.1 x 10 43 


4.9 x 10 42 


3.0 x 10 43 



Notes. Summary of the SED parameters for the two Fenni/LAT spectra 
considered (low-energy threshold equal to 300 MeV and 1 GeV) and 
for the two cases proposed in the Swifr/XRT analysis (case A and case 
B). For all the cases considered, the minimum and maximum Lorentz 
factors of the electron distribution are set to y e ^ m = 300 and •y (Mna x = 
5 x 10 6 respectively; the angle to the line of sight is 6 = 1°, the Doppler 
factor is 6 = 30, and the index of the electron distribution after the break 
is a 2 = 4.0. The electron-distribution normalization parameter K is in 
units of cm" 3 ; the magnetic field B is in G; the emitting-region size R is 
in cm; the energy densities u e B are in erg cm" 3 ; the variability timescale 
T var is in hours and the jet luminosity Lj e t is in ergs s 



4. SSC modelling of the SED 

The non-simultaneous SED of 1RXS J101015.9-31 1909, cor- 
rected for Galactic absorption, is shown in Fig. [7] Historical data 
taken from NEEO are also shown. Before 2006 the source has 
been observed in radio, infrared, optical and X-rays. A discus- 
sion of the accuracy of the cross-calibration between Fermi/LAT 
and H.E. S.S. (evaluated at 4% based on the Crab nebula) can be 
found in lMeveretaTI (l2010h . 

The optical flux contribution from the host galaxy has been 
evalua ted using data from the 2MASS Extended Source cata- 
logue dJarrett et alJfe OOO). Based on the magnitude values eval- 
uated for different radii of aperture (r > 5") and the effective 
radius of the galaxy (r e jf = 3.03" in the J band), we estimate 
this contribution in a 4" radius of ap erture a s m sa \ 14.3 for 
the magnitude in the J band, following lYound (119761) . The mag- 
nitude obtained has been used to properly rescale the template 
of a giant elliptical galaxy spectr um (evaluated using PEGASE) 
dFioc & Rocca-Volmerangell997l) . 

As shown in Fig. [7] in infrared light, the host galaxy domi- 
nates the AGN emis sion. This is cons i stent w ith the optical spec- 
trum measured bv lPiranom onte et aT1 (l2007l) when evaluating the 
redshift of the source (see Fig. 2 and A. 1 in their paper) and 
with the fact that the variability amplitude in the B band is sig- 
nificantly larger than that in the R band. 

The emission from the ac tive nucleus is de scribed using a 
stationary one-zone SSC code dKatarzvhski et al.l2001l) : a spher- 
ical plasma blob (characterised by its radius R) moving with 
Doppler factor 5 in the relativistic jet (with being the angle 
to the line of sight) is filled with a homogeneous magnetic field 
B and a stationary, non-thermal electron distribution. The syn- 
chrotron emission from these electrons is responsible for the low 
energy bump, peaking in the X-ray band, and is then Compton- 
upscattered by the electrons themselves, to produce the y-ray 



11 NASA/IPAC extragalactic database, http://ned.ipac.caltech.edu 
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emission. Pair production (y + y — > e~ + e + ) inside the blob 
is not negligible, and is taken into accoun t using the cross sec- 
tion evaluated by lAharonian et al.l (120081) . The interaction be- 
tween VHE photons escaping the source and the infrared extra- 
galactic background light (EBL) produces e + e~ pairs as well, 
and induces an absorption in the observed VHE spectrum of the 
source. Thi s effect has been taken in to account using the EBL 
model from lFranceschini et al.l (120081) . which is compatible with 
the EB L limit inferred from VHE observations (Aharonia n et al.l 
2006a). The primary electron distribution (N(y e ), where y e is 
the Lorentz factor of the electrons), defined between y e m i„ and 
Tf.max, is modeled using a broken power-law function with 
normalization K (defined as the number density of electrons at 
y e — 1, in units of cirT 3 ) and indices a\ below, and above the 
break Lorentz factor 7V,b re ak- 

The two previously-mentioned X-ray spectral hypotheses 
(assuming an intrinsic break or an additional absorption) have 
been considered as lower and upper limits for the synchrotron 
emission from the blob. In case A, the synchrotron peak energy 
corresponds to the observed X-ray break energy, while in case B, 
the synchrotron peak falls between UV and X-rays. Whereas in 
case A the emission from the blob cannot explain both the X-ray 
and optic al/UV data, in case B the synchrotron component, to- 
gether with the emission from the host galaxy, can reproduce the 
infrared to X-ray observations. In both cases, the historical radio 
data are not taken into account, as it is more likely that they are 
produced in the extended jet. In order to study how the evalua- 
tion of the Fermi/LAT slope (which depends on the low-energy 
threshold, as described in Sec. 3.1) affects the overall SED, the 
modelling has been performed for the two Fermi/LAT spectra 
evaluated above 300 MeV and 1 GeV For simplicity, only the 
modelling of the SED with the Fermi/LAT spectrum evaluated 
above 1 GeV is presented in Fig. [7] 

The minimum and maximum Lorentz factors of the electron 
distribution cannot be constrained by the data, and they have 
been fixed at y esr ( m = 300 and y e , mdx = 5 x 10 6 . The index of the 
electron distribution after the break a-i is completely constrained 
by the observed X-ray photon index above the break, and has 
been fixed at 4.0. The slope a\ is constrained by the Fermi/LAT 
photon index (for the two cases, A and B) and by the optical/UV 
data points (only for case B), and it has been fixed at 2.2 and 2.0 
for a Fermi/LAT spectrum evaluated above 300 MeV and 1 GeV, 
respectively. 

A good description of the SED can be obtained assuming an 
angle to the line of sight 6 — 1° and a Doppler factor 6 = 30, cor- 
responding to a bulk Lorentz factor of 16. The other free parame- 
ters (B, R, 7e t break, K) are different between the cases considered, 
and their values are indicated in Tab. [5] together with the evalu- 
ation of the electron energy density u e = mc 2 J dy e y e N(y e ) and 
the magnetic energy density in the blob ub = B 2 l%n. The u e /us 
value is higher in case A (intrinsic break, and higher synchrotron 
peak energy) than in case B (additional absorption effect, and 
lower synchrotron peak energy), reflecting the fact that the ratio 
between the inverse Compton and the synchrotron component is 
higher in the first case. The lower limit on the variability time- 
scale, evaluated for the emitting region size and the Doppler fac- 
tor assumed in the modelling, roughly corresponds to 1 and 25 
hours for the cases A and B, respectively, consistent with the 
variability time-scale observed by ATOM. 

The difference between the Fermi/LAT spectra evaluated for 
different energy thresholds affects the evaluation of the electron- 



12 For a j ustific ation of the use of this hypothesis see, for example, 
lKirketalJ ( fT998l) . 



distribution slope before the break (ai), and, consequently, it in- 
duces a variation on the normalization factor K and on the break 
Lorentz factor y ei break, modifying the value of the electron energy 
density inside the emitting region. 

In case A, the observed flux at low frequency (infrared to 
UV) cannot be explained by the emission of the blob. An ad- 
ditional component is required (not shown in FigJT}, such as the 
emission from the extended jet, dominating the non-thermal con- 
tinuum from radio to UV, and being responsible for the variabil- 
ity observed in ATOM data. 

On the other hand, in case B, the low-frequency emission can 
be described by the blob-in-jet component plus the contribution 
from the host galaxy, with the former being at the origin of the 
observed optical variability. 

It should be noted that, in this case, the UV flux is slightly 
underestimated. To better describe the data, the model would 
need a harder slope a\ in apparent conflict with the GeV con- 
straints derived from Fermi/LAT. However, the uncertainties 
(both statistical and systematic) on the evaluation of the GeV 
slope, as well as on the value of the second absorber in case B 
(the error on the Eb-v value used for the dereddening of the data 
is about 20%; this uncertainty has not been taken into account in 
the plotting of the SED), can still explain this discrepancy. 

As mentioned above, the two cases discussed here (A and 
B) are best considered as lower and upper limits for the SSC 
blazar emission. The real scenario may be more complex and lie 
between these two limiting cases. 



5. Conclusions 

The blazar 1RXS J101015.9-31 1909 has been observed by 
H.E.S.S. between 2006 and 2010, leading to the discovery of 
its VHE emission with a significance of 7.1 standard deviations. 
The time-averaged VHE spectrum of this blazar is soft, with a 
photon index of T = 3.08 + 0.42 stat ± 0.20 sys and a flux 0.8% of 
that of the Crab nebula. The detection has been made using more 
powerful analysis methods, which provide an enhanced sensitiv- 
ity at lower energies. Observations at other wavelengths have 
been analysed in order to have a multi-band view of the SED 
of this AGN newly-detected in the VHE range. In particular, a 
careful analysis of the HE emission in Fermi/LAT data reveals a 
detection of this AGN at a significance level of about 8.2 stan- 
dard deviations. 

No detectable variability is seen in the available VHE 
(H.E.S.S.), HE (Fermi/LAT) or X-ray (Swift/XRT) datasets, the 
latter being the most constraining for the evaluation of the vari- 
ability at the * 10 3 s timescale of the X-ray observations that 
have been analysed. On the other hand, the optical data from 
ATOM, contemporaneous to the H.E.S.S. observations, vary in 
both the B and R filters. The fitted position of the VHE emission 
is consistent with the positions of the X-ray and radio sources. 

The study of the SED implies that 1RXS J101015.9-31 1909 
is an HBL with the usual double-peaked SED, where the high- 
energy bump is fully characterized thanks to the detection of the 
source in the H.E.S.S. and Fermi/LAT datasets. The physical 
origin of the spectral break observed in the X-ray spectrum can- 
not be firmly determined, and two hypotheses, one assuming an 
intrinsic break and another considering an additional absorption 
effect, have been presented. The corresponding X-ray spectra are 
considered as the lower and upper limits of the intrinsic blazar 
emission at that wavelength. The SED can be reproduced rea- 
sonably well using a stationary one-zone SSC model, with pa- 
rameter values compatible with those commonly assumed for 
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Energy (eV) 




v(Hz) 

Fig. 7. Spectral energy distribution of 1RXS J101015.9-31 1909. The H.E.S.S. spectrum is represented by the green bow-tie at 
the highest energies. The measured Fermi/LAT spectrum above 1 GeV is represented by the orange full bow-tie, while the binned 
spectral points or upper limits are shown with empty circles (see Sec. 13. II for details). The Swift/XRT spectrum in case A is shown 
with green diamonds, while case B is represented with pink crosses. NED database archival data are shown with grey asterisks, 
while the grey cross represents the flux measured by the SUMSS radio survey (see Sec. [TJ. Swift/UVOT data are shown with 
light-red stars; ATOM data are shown with light-blue squares. Deabsorbed infrared-to-UV data are plotted for case B only to avoid 
cluttering. The total SSC emission model (including absorption by the EBL) in case A is represented by the solid line, while the 
modelling in case B is represented by the dashed line. The intrinsic emission models (corrected for EBL absorption) are presented 
by the thinner lines at higher VHE flux. A template of a spectrum of a giant elliptical galaxy is also shown in the plot by the dashed 
purple line in the optical range. 



relativistic jets. However, if the origin of the X-ray break is in- 
trinsic, the model requires another component (most likely the 
emission from the extended jet) in order to describe optical and 
UV data. 

Current uncertainties on the intrinsic X-ray spectrum, to- 
gether with the non-simultaneity of the infrared to UV data, re- 
main limiting factors for a realistic modelling of the low-energy 
bump. Future simultaneous multi-wavelength observations are 
required in order to determine whether optical/UV and X-ray 
photons are produced in the same emitting region, thus con- 
straining the origin of the observed X-ray break. 

The high-energy part of the SED turns out to be better con- 
strained, with modelling only limited by the current uncertainty 
in the Fermi/LAT spectral index. Given the surveying strategy 
of Fermi/LAT, further data are expected to allow a more precise 
modelling of the source to be attained. 
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